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Abstract 

A recent HyperCP observation of three events in the decay S + — ► u~ is suggestive of a new particle 
with mass 214.3 MeV. In order to confront models that contain a light Higgs boson with this observation, 
it is necessary to know the Higgs production rate in hyperon decay. The contribution to this rate from 
penguin-like two-quark operators has been considered before and found to be too large. We point out that 
there are additional four-quark contributions to this rate that could be comparable in size to the two-quark 
contributions, and that could bring the total rate to the observed level in some models. To this effect we 
implement the low-energy theorems that dictate the couplings of light Higgs bosons to hyperons at leading 
order in chiral perturbation theory. We consider the cases of scalar and pseudoscalar Higgs bosons in 
the standard model and in its two-Higgs-doublet extensions to illustrate the challenges posed by existing 
experimental constraints and suggest possible avenues for models to satisfy them. 
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I. INTRODUCTION 



Three events for the decay mode S + — > p/x + /i~ with a dimuon invariant mass of 214.3 ±0.5 MeV 
have been recently observed by the HyperCP Collaboration [1]. It is possible to account for these 
events within the standard model (SM) when long-distance contributions are properly included [2]. 
However, the probability of having all three events at the same dimuon mass in the SM is less than 
one percent. This suggests a new-particle interpretation for the events, for which the branching 
ratio is (3.1±?;| =b 1.5) x 1(T 8 [1]. 

This possibility has been explored to some extent in the literature, where it has been shown that 
kaon decays place severe constraints on the couplings of the hypothetical new particle [3-5]. In 
particular, it was found that the flavor- changing coupling of the new state, X, to ds has to be of 
a pseudoscalar or axial- vector nature to explain why the state has not been seen in K — > 7r/i + /x~. 
At least one model containing a particle with these properties has appeared in the literature [6]. 

All these previous analyses of X considered only the effects of two-quark operators for dsX. 
However, it is well known in the case of light-Higgs production in kaon decay that there are also 
four-quark operators that can contribute at the same level as the two-quark ones [7-10]. These 
four-quark contributions are most conveniently described in chiral perturbation theory (%PT) which 
implements low-energy theorems governing the couplings of light (pseudo)scalars to hadrons. In this 
paper we generalize existing studies appropriate for kaon decay to the case of hyperon decay. This 
allows us to discuss the production of light (pseudo)scalars in hyperon decay consistently, including 
the effects of both the two- and four-quark operators with the aid of xPT. We consider the cases 
of scalar and pseudoscalar Higgs bosons in the SM and in the two-Higgs-doublet model (2HDM), 
expressing our results in a form that can be easily applied to more complicated Higgs models. 

This paper is organized as follows. We begin by collecting in Sec. II the existing constraints on 
light Higgs bosons from kaon, S-meson, and hyperon decays if we interpret the HyperCP events 
as being mediated by a light Higgs boson. In Sees. Ill and IV we compute the production rates in 
both kaon and hyperon decays for a light scalar and pseudoscalar Higgs boson, respectively. Finally 
in Sec. V we summarize our results and state our conclusions. 



II. SUMMARY OF EXISTING CONSTRAINTS 



In Ref. [3] we parameterized the possible couplings of the new particle, X, to ds and /x/x assuming 
that it had definite parity. Whereas this is a reasonable assumption for the diagonal couplings of 
X to fermions, it is not for its flavor-changing neutral couplings (FCNCs). Two-quark FCNCs are 
predominantly induced by Higgs-penguin diagrams, which result in left- and right-handed couplings, 
implying that the scalar and pseudoscalar ones are present simultaneously. For this reason, we revisit 
the existing constraints for X being a scalar particle, H, or a pseudoscalar particle, A, assuming 
them to have two-fermion FCNCs described by 

C Hsd = 9 ~f K + 7 5 > + m d d(l - 75 )s] H + H.c. , (la) 
J^Asd = — [™s d(l + l 5 )s - m d d(l - 75 )s] A + H.c. , (lb) 
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where the g's are coupling constants, m q is a quark mass, and v = 2" 1 / 4 G F = 246 GeV. In 
addition, the diagonal couplings to charged leptons are assumed to have definite parity and be 
proportional to the lepton mass, 

£ M = 9 -^MH, C M = l -^hdA. (2) 

For a (pseudo) scalar of mass 214.3 MeV, it is then natural to assume that the decay X — > 
will dominate over the other kinematically allowed modes: X — > e + e~, vv, 77. We will restrict 
ourselves to this case, assuming that B(X — > ~ 1. This is true, for example, for a light 

SM Higgs boson where g e — 1, or for light pseudoscalars in the 2HDM types I and II, where 
g e = cot (5 and — tan /3, respectively. In all these cases decays, X — > e + e~ are suppressed at least 
by {mjm^f ~ 10" 5 . 

To be consistent with the HyperCP observation, X must be short-lived and decay inside the 
detector. This is compatible with the estimate for the total width ~ 10 _7 MeV [5] of a pseu- 
doscalar particle, A. It was shown in Ref. [3] that the muon anomalous magnetic moment imposes 
the constraint 

\g e \ < 1.2 . (3) 
A coupling satisfying this constraint implies a width 

1^ < 3.7 x 1(T 7 MeV , (4) 

consistent with the observation. In contrast, the corresponding constraint for a scalar particle is 
\g e \ < 0.98, leading to a longer lifetime, 

T H < 6.9 x 1(T 9 MeV . (5) 

The estimated lifetime for the HyperCP particle is therefore consistent with that of a pseudoscalar 
or scalar that decays predominantly into muons. 

In addition to the two-quark contributions to the amplitudes for K — > ttTI(A) and S + — > pTi(A) 
induced by the interactions in Eq. (1), we will also include contributions arising from the usual SM 
four-quark \AS\ = 1 operators, along with flavor- conserving couplings of 7i{A). We will adopt the 
chiral-Lagrangian approach to evaluate the hadron-level interactions. 

Later on we will discuss specific models and consider the bounds appropriate for them, including 
all the relevant two- and four-quark contributions. It is useful to start with one example to illustrate 
the ingredients needed to construct a model that can satisfy all the existing constraints. For this 
purpose, we consider a pseudoscalar A with two-quark couplings as in Eq. (lb) supplemented with 
simple parameterizations for the four-quark amplitudes for both kaon and hyperon decays. For 
£>-meson decay, we assume that the two-quark contribution completely dominates. 
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A. K — ► 7vA 



Introducing the dimensionless quantity M±k for the four-quark contribution, we express the 
amplitude for K ± — > n^A and its branching ratio, respectively, as 

iM(K± - ttU) = g A ^< - M 4X ^ , 

v v 

B{K± -> tt ± ^) = 4.43 x 10 s - 1.08 M 4K \ 2 . (6) 

This mode is constrained by its nonobservation in the BNL E865 [11] or FNAL HyperCP [12] 
measurements of K ± — > n ± fi + fi^. It is also constrained by its nonobservation in the isospin- 
related mode K$ — > ir°fi + fi~ by CERN NA48 [13]. Of these three experiments, E865 had the best 
statistics, collecting 430 events in K + — > 7r + /i + /i~. A new particle A of mass 214.3 MeV would 
have contributed only in their first dimuon-mass bin, where 0.21 GeV < < 0.224 GeV and 
approximately 30 events were observed. To obtain a conservative bound, we assume that all the 
events in the first bin are statistically Gaussian and can be attributed to the new particle (either a 
scalar or a pseudoscalar) . Further assuming uniform acceptance, we obtain at 95% C.L. 

B(K + -> 7T+X) < 8.7 x 10~ 9 . (7) 

The NA48 Collaboration collected 6 events for K s — > 7r /! + /i~ [13], and none of them have the 
214.3-MeV invariant mass required if they originate from the new particle A. Using the Ks flux 
and the acceptance at low m in Ref. [13], we estimate a single event sensitivity of (5.3^o 4) x 10~ 10 . 
With no events observed and Poisson statistics, this translates into the 95%-C.L. bound 

B(K S -> 7r°X) < 1.8 x 10~ 9 . (8) 

We employ these bounds when we discuss specific models, but for now we use the E865 result in 
Eq. (7), combined with Eq. (6), to find 

|^-1.08M 4X | < 4.4xl0~ 9 . (9) 



B. £+ -> pA 

In this case, we need two new dimensionless quantities A A and _B 4 to parameterize the effect of 
the four-quark operators, writing the amplitude as 

- pA) = ip (A pA - B pAlb ) S+ , (10a) 

where 

_ m^-m N f 

1 V V 

B A = g A (D-F) m * + mN f K 2 +B 4 ^, (10b) 
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the parameters D and F coming from a chiral Lagrangian to be discussed in a later section and 
/„. = 92.4 MeV being the pion-decay constant. The resulting branching ratio is 

B(Z + ^pA) = 1.91 x 10 6 1^ + 0.36 A 4 \ 2 + 4.84 x 10 4 \g A + 0.U B 4 \ 2 (11) 

with the choice D — F = 0.25. Combining the statistical and systematic errors of the HyperCP 
measurement [1] in quadrature, we require 

-> pA) = (3.1111) x 1(T 8 , (12) 

and therefore 

\g A + 0.36 A 4 \ = (1.3 ±0.6) x 10~ 7 , (13) 

where we have used the larger of the errors in Eq. (12) and ignored the contribution from the 
P-wave term in Eq. (11), assuming that B4 < A4. This assumption is satisfied by all the models 
we discuss, but when checking a specific model, we do so without neglecting £> 4 . 

A comparison of Eqs. (9) and (13) shows why it is not possible to have a (pseudo) scalar with 
penguin-like flavor- changing neutral couplings, as in Eq. (1), as an explanation for the HyperCP 
result given the constraints from kaon decay. It also shows how this is no longer true if there are 
four-quark contributions to the amplitudes that are comparable to the penguin amplitudes. In 
particular, if we assume that in a given model g A , M 4K , and A 4 have comparable magnitudes, we 
see that in order to satisfy both Eqs. (9) and (13) we need a cancelation between the two- and 
four-quark contributions to the kaon amplitude that reduces them by a factor of about 20. As we 
will show in later sections, this is possible in many models. For this cancelation to work, however, 
g A and M 4 k must also have similar phases. As we will see, this is a requirement that is much harder 
to satisfy. In the simple models we consider in this paper, the phase of g A is much larger than the 
phase of M^k so that the cancelation does not happen for the imaginary part. 



C. b^sX 

Finally we consider the constraints on the new particle from its nonobservation in £?-meson 
decay. In this case, the four-quark contributions are negligible, and we can neglect m s compared to 
m b . The Lagrangian for b — > sX can then be expressed as 

C Xbs = 9 -^s{l+ l5 )bX + H.c. , (14) 
where g' = g' n (ig' A ) for X = H (A). This leads to the partial decay rate 

r(»-rf) = lsf|i. (15) 

Using for illustration m b = 4.3 GeV and the B + lifetime [14] results in 

B{b^sX) = 1.3 x 10 8 |#'| 2 . (16) 
One could obtain a similar number for b — > dX. 
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The latest experimental average B(b — > sfi + fi~) = (4.27l};22) x 10~ 6 [14] covers the full kinematic 
range for m . To constrain g', it is better to limit the comparison to the measured rate at the 
lowest measured invariant-mass bin. BABAR quotes in Table II of Ref. [15] 

B{b -> ^O m ^_ e[0 . 2GcViL0GeV] = (0.08 ± 0.3618;8I) x 10- 6 . (17) 

This is an average for electrons and muons, but no noticeable difference between them was found. 
Belle quotes on Table IV of Ref. [16] the corresponding number 

B(b -> s£ + n me+e _ e[0 , 2GcVA , 0GcV] = (11.3 ± 4.81*;?) x lO" 7 . (18) 

To be conservative, we constrain the Higgs coupling by requiring that the induced rate be below 
the 95%-C.L. upper-range of the measured b — > s£ + £~ rate in the lowest measured bin. Thus, 
combining errors in quadrature for the more restrictive BABAR result gives 

B (b -> sl+r) me+e _ <1GcY < 8.0 x 10- 7 (BABAR) (19) 

and correspondingly 

\g'\ < 7.8 x 10~ 8 . (20) 

The exclusive B — > (K, K*)fi + fi~ modes have been measured, but the resulting constraints are not 
better than Eq. (20). This constraint, Eq. (20), is difficult to satisfy in models where g A and g' are 
related by top-quark CKM angles, as happens in the simple models we consider here. 

III. SCALAR HIGGS BOSON 

In this section we discuss in detail the case of a light Higgs boson in the standard model and in 
the two-Higgs-doublet model. We will use known low-energy theorems to implement the four-quark 
contributions to kaon and hyperon amplitudes. 

A. Two-quark |AS| = 1 interactions 

The effective Lagrangian for the sdH coupling, where Ti is either the standard-model Higgs 
boson H° or the lightest scalar Higgs-boson h° in the 2HDM, has been much discussed in the 
literature [7, 8, 17-19] and can be written as Cnsd in Eq- (la), where 

9n = yjfe E nXtV^Fiq), (21) 

q=u,c,t 

with V kl being the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and F(q) depending 
on the model. In the SM, for a Higgs mass much smaller than the W mass, 

F(q) = 3/4 , (22) 

whereas in the 2HDM the expression for F(q) is much lengthier [18, 19]. 
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Using CKM and mass parameters from Ref. [20], we find in the SM 

9n = (-1.3 - 0.6i) x 1(T 6 , (23) 

to be compared with Eqs. (9) and (13) above. Employing the expression for F(q) derived in Ref. [19], 
we obtain a similar number in the 2HDM type II, for instance, 

9n = (5.0 + 1.9i) x 10~ 7 (24) 

for the parameters 

tan/3 ~ 2.57, sin(/3 - a) ~ 0.149, m H+ = 250 GeV , (25) 

where tan j3 is the ratio of vacuum expectation values of the two Higgs doublets, a the mixing angle 
in the neutral-Higgs-boson mass matrix, and mu+ the mass of the charged Higgs bosons. 1 We note 
that the a and (3 values above satisfy the constraint sin 2 (/3 — a) < 0.06 from LEP [21]. We see 
right away that g H can be in the right ball park to explain the HyperCP observation, Eq. (13), but 
conflicts with the kaon bound, Eq. (9). 

To evaluate the hadronic amplitudes from this 2-quark contribution, we employ chiral perturba- 
tion theory. Using the operator matching of Ref. [3] , we write the lowest-order chiral realization of 
C Hsd as 

C H = b D (B{h n ,B}) + b F (B[h n ,B])+b (h n )(BB) + \f 2 B,{h H ) + H.c. , (26) 
where (•••) = Tr(- • • ) in flavor-SU(3) space, f = f n = 92.4 MeV, and 

h H = -2g n (tfhM? + tMh£) ™ , (27) 

with h being a 3x3-matrix having elements h kl = S k2 S 3l which selects out s — > d transitions, 
M = diag(m, to, m s ) = diag(ro 2 , to 2 , 2m\ — m%)/(2B ) the quark-mass matrix in the isospin- 
symmetric limit m u = m d = to, and the baryon and meson fields represented by the usual 3x3- 
matrices B and E = = e l(f ^ , respectively. 

To derive amplitudes, we also need the chiral Lagrangian for the strong interactions of the 
hadrons [22, 23]. At leading order in the derivative and m s expansions, it can be written as 

4 = (Bt^(d^B+[V^B]))-m (BB) + D(B^ l5 {A^B}) + F(B^ l5 [A^B]) 
+ b D (B {M + , B}) + b F (B [M+, B])+b (M+> (BB) 

+ \f <^Et ^E) + \fB, (M + > , (28) 

where V M = \ (t > + <9 M £) , m is the baryon mass in the chiral limit, A^ = \ (£ d^tf - <9^£) , 
and M + = ^M^ + £M f £, with further details being given in Ref. [3]. 

From C-u and C s , we derive the leading-order diagrams shown in Fig. 1 for E + — > pTC, yielding 
the amplitude 

H " mi- — m z 



K 7T 



2 ™2 



- 9n( D ~ F ) 2 T P^5^ , 29 

v m A K — 



1 We have also set k = ■m 2 H+ /m% v in F(q), where k is defined in Ref. [19]. 
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(a) 



FIG. 1: Diagrams contributing to S 



(b) 



pTL arising from C-usd at leading order in xPT. The square 



vertices come from L-yi in Eq. (26), and the solid vertex from C s in Eq. (28). 



where the two terms correspond to the two diagrams, respectively, m s N are isospin-symmetric 
masses, and we have used the relations m E — m N = 2 (b D — b F ) (m s — m) , m 2 K = B (m + m s ) , and 
m 2 = 2B m derived from C s . Numerically, we will allow D and F to have the ranges 0.6 < D < 0.8 
and 0.4 < F < 0.5 [23], leading to 

0.1 < D — F < 0.4 , (30) 

which is their combination occurring in our amplitudes. 

It follows that the contribution of C-H S d to the branching ratio of S + — > pTi for m n = 214.3 MeV 
and the middle value D — F = 0.25 is in the SM 

B 2q (J: + -> pH) = (40 + 1) x 10" 7 , (31) 

where we have ignored the imaginary (CP violating) part of the amplitude, and the two num- 
bers correspond to the contributions from the scalar and pseudoscalar flavor- changing couplings, 
respectively. Evidently, the scalar contribution is much larger than what HyperCP saw, but the 
pseudoscalar contribution is within the range. This, however, is only part of the story, as there are 
in addition 4-quark contributions to be discussed in the next subsection. 

Also from we derive the leading-order diagram for K — > irTi, which is that in Fig. 1(a) with 
S + (p) replaced by K (tc) and arises from the scalar coupling in Cnsd- The resulting amplitude is 

M 2q (K + - n + H) = -V2M 2q (K° - 7r°H) = ~ g ^ K , (32) 

and so M.2q (K L — > tt H) = — ReA4 2q (K + — > n + H). Dropping again the imaginary parts of the 
amplitudes, we obtain in the SM the branching ratios 

B 2q (K + -> tv + H) = 9.3 x 10^ 4 , B 2q (K L -> tv°H) = 3.9 x 10~ 3 . (33) 

These numbers would easily be incompatible with that in Eq. (7) and the 95%-C.L. bound 2 

B(K L -> ttVV - ) < 4 -9 x 10" 10 , (34) 

but, as in the S + case, there are 4-quark contributions that have to be considered as well. 



2 We have inferred this number from Ref. [24] which reported B(Kl — > n°fi + fi ) < 3.8 x 10 10 at 90% C.L. 
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The situation is similar in the 2HDM. Adopting the real part of the coupling in Eq. (24), for 
example, we find 

B 2q {^ + -> pH) = (56 + 1) x 1(T 8 , 
&2q( K+ -" = 1-3 x 10 -4 , B 2q (K L -> 7T ^) = 5.4 x 10 -4 . (35) 



B. Four-quark |A5| = 1 interactions 



The hadronic interactions of a light Higgs-boson due to 4-quark \AS\ = 1 operators are best 
accounted for in the chiral Lagrangian approach. The dominant contribution is generated by the 
| A/ 1 = \ component of the effective Hamiltonian transforming as (8 L , 1 R ). The corresponding 
Lagrangian at leading order is given by [23, 25]. 



h D (B{?hS,B}) + h F (B[?h£,B]) 

+ V 2 (hd^Ptf) + 2%f 2 B (hZM+tf) + H.c. 



(36) 



where h D F can be extracted from hyperon nonleptonic decays, 7 8 = —7.8 x 1CT 8 from K — > nn, 
the sign following from various predictions [7-9, 26], and 7 g is unknown as it does not contribute 
to any process with only kaons and pions. 

The 4-quark |A5'| = 1 interactions of a light Higgs-boson TC arises from its tree-level couplings 
to quarks and bosons, as well as from its coupling to gluons induced by a triangle diagram with 
heavy quarks in the loop. To obtain the relevant chiral Lagrangians, one starts with £ SjW above and 
follows the prescription given in Refs. [7-9, 19]. The results are 



= (J c, f (9>tf + i C2 f Bo {M+ ) + I fB (M + - M + >) ^ - k ± m (BB) 



n 



+ 



k 2 (b D (B {M + ,B}) + b F (B [M + , B]) + b (M + ) {BB)) 



n 

V 



(37) 



C n = 



where 



with 



7 8 c 3 f (hd^ + 27 8 c 4 f 2 B (hZM+?) + 27 8 f 2 B (h£(M + - M + )tf) 



n 



J V 



n 



+ k 3 (h D (B{?ht,B}) + h F (B[?ht,B]))- + H.c, 



(38) 



2k@ , c 2 — + 1 , c 3 — 2/c^y , c 4 — 5/^ 2/c^ + 1 , 

^i = ' ^2 = 1 ) ^3 = 3&G- — 2/c w , (39) 



^g = 2(2/ ^+ ^ , M = diag(/c>, fc d m, fc d m s ) , 



(40) 
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the expression for k G corresponding to 3 heavy and 3 light quarks. The parameters k ud , k w , and 
k G come from the couplings of 7i to light quarks, W ± , and the gluons, respectively, and depend on 
the model of the Higgs sector. Thus 



A , 



k u — kd — 



cos a 



k u — kd 
cos a 



k 



w 



in the SM , 



sin P ' 
k A = - 



k w = sin(/3 — a) 



in the 2HDM I , 



sin a 



k 



w 



sin(/3 — a) 



in the 2HDM II 



(41) 
(42) 

(43) 



sin/5 ' a cos/5 ' 

The parameters c x 2 3 4 for the meson terms have already been obtained in the literature [7-9, 19, 27], 
whereas the new ones k 123 follow from how the baryon parameters depend on masses: m ~ A, 
b D F0 ~ 1, x+ ~ Am g , and h DF ~ A 3 /m^/, where A is a QCD mass scale. Note that we work in 
that basis in which the mass terms in the Lagrangians are not diagonal and must therefore include 
the corresponding tadpole diagrams in our calculation. 

For S + — > pH, we derive from the diagrams shown in Fig. 2, finding 



M 4q (Z + ^pH) = (k d -3k G + 2k w ) 



h D -h, 



pE+ 



+ 4(k G -k w )(D-F)^^±^L 



in 



K 



m: 



m 2 K -m 2 n 



(44) 



where the first term comes from the upper three diagrams, which are at leading order, and the 
second term results from the lower two diagrams, which are at next-to-leading order. Now, the 
combination ho— hp also occurs in the amplitude for S + — > pn , which we write as 



where from £<, 



M(Z + ^pn°) = ip (A pn0 - 5 p7r075 ) £+ , 
-h D + h F h D -h F m^ + m N 



2/ 



(D-F)- 



2/ 



m s -m N 



Since from experiment [28] 



-3.25 x 10" 



26.67 x 10" 



(45) 

(46) 

(47) 
pH we will 



up to an overall sign, in our numerical evaluation of the 4-quark contributions to E + 
explore different h D —h F values accordingly. 

We can also derive from Cl,J the corresponding leading-order diagrams for K — > ixTi, which are 
the upper three in Fig. 2 with S + (p) replaced by K (jr) and yield 

7s 



M 4q {K+^n + H) = f[2(k w -k G ){m 2 K + ml-m 2 n ) + (k d -k u )ml] 



+ ^4(k G 



k w ) m 2 K , 



(48) 



M 4q (K° -> 7i°n) 



7s 

% 



m? m 



K 



m 



K 



mi 



+ 



V2i 



4(k w - k G ) m\ + (k d - k u ) 



mt m 



K 



(49) 
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FIG. 2: Diagrams contributing to S + — > arising from the 4-quark operators. The square vertices come 
from £^ in Eqs. (36) and (38), whereas the dots are from £g in Eqs. (28) and (37). 



Since 7 8 is unknown, in evaluating its effect on K + — > 7r + ?^ we will allow it to vary from —10 to 
10 times 7 8 . Naively we would expect 7 8 and 7 8 to be of the same order. 



C. Total contributions 



The total amplitude for K + — > -K + 7i comes from the sum of the contributions in Eqs. (32) 
and (48). If the CP- violating terms in the amplitudes are ignored, it is possible for the 2-quark 
and 4-quark contributions to cancel. We show this possibility in Fig. 3, where we plot the resulting 
branching ratio in the SM as a function of the ratio r 8 = 7 8 /7 8 for m n = 214.3 MeV. We find 
that B(K + — > 7i + 7i) = when r 8 ~ —5.1 and that, as the figure indicates, for only a very narrow 
range of r 8 around this value does the branching ratio ever fall below the upper limit in Eq. (7). In 
Fig. 3, we also plot the corresponding branching ratio of the isospin-related mode Kl — > 7r°7i. 

For S + — > pTi, the total amplitude results from adding the contributions in Eqs. (29) and (44). 
Including only the real part of amplitudes again, and using r 8 ~ —5.1 determined above, we plot 
in Fig. 4 the branching ratio in the SM as a function of D—F for the range in Eq. (30). This figure 
shows that the curve resulting from the P-wave fit using Eqs. (46) and (47) satisfies the HyperCP 
constraints for certain D—F values. 

In Figs. 5 and 6, we display the corresponding branching ratios in the 2HDM II obtained using 
the parameters in Eqs. (24) and (25). In contrast to the SM case, here B(K + — > TT + 7i) = when 
r 8 ~ 6.7, but the vanishing of the Kl rate occurs at a different r 8 value due to Ai^Ki — > 7i°H) 
and — ReA4 4:q (K + — > 7r + H) being unequal with k u ^ k d in Eq. (43). 3 As a consequence, the two 
kaon constraints cannot be satisfied simultaneously. Furthermore, the S + — > pH curve that falls 
within the HyperCP limits is the one resulting from the S-wave fit using Eqs. (46) and (47). 

To summarize this section, a light Higgs boson in the SM can be made compatible with the 
empirical bounds for S + — > pTi, while satisfying the constraints from K — > irli., if the real part 



3 We note that, although 7 8 is not known from experiment, there are model calculations [8, 26] of it yielding 
|7 8 /7 8 | ~ 0.2 . This would make the kaon rates greatly exceed their bounds, as can be seen from Figs. 3 and 5. 
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FIG. 3: Contributions of real parts of total amplitudes for K + — * tt + 7{ (solid curve 
(dashed curve) in the SM to their branching ratios as functions of r 8 = 7g/7g for 
horizontal lines are the corresponding upper bounds in Eqs. (7) and (34). 



and K L -> ir°H 
-- 214.3 MeV. The 
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FIG. 4: Contribution of real part of total amplitude for S + — ► pTC to its branching ratio in the SM as 
function of D—F for = 214.3 MeV and r g ~ —5.1. The solid (dotted) curve corresponds to hj^—hp 
extracted from the P-wave (S-wave) fit to the S + — > pn data using Eqs. (46) and (47). The dashed lines 
correspond to the upper and lower bounds in the Hyper CP result. 



of the 2-quark (penguin) contribution to the respective amplitudes is combined with the 4-quark 
contribution. Moreover, in the 2HDM such a particle can satisfy all these constraints if its diagonal 
couplings to the up- and down-type quarks are the same. For this to happen in either model, it 
is necessary for the two amplitudes to cancel precisely, and we have shown that this is possible for 
certain values of the hadronic constants 7 8 , hp— hp, and D—F . Although 7 8 is not known, unlike 
hp,— hp and D—F which are extractable from hyperon nonleptonic and semileptonic decays [23], it 
has a definite value in the SM and cannot be fine-tuned. We note that in all the S + — > pFC cases 
discussed above the P7 5 E + term in the amplitude is small compared to the pT, + term and that, 
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FIG. 5: Contributions of real parts of total amplitudes for K + — » tt + TL (solid curve) and Kl — ► 7T°Ti 
(dashed curve) in the 2HDM II to their branching ratios as functions of r 8 = 7§/7g for = 214.3 MeV 
and the parameters in Eq. (25). The horizontal lines indicate the upper bounds in Eqs. (7) and (34). 
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FIG. 6: Contribution of real part of total amplitude for S + — ► p7i to its branching ratio in the 2HDM 
as function of D—F for = 214.3 MeV and the parameters in Eq. (25). The solid (dotted) curve 
corresponds to hp— hp extracted from the P-wave (S-wave) fit to the S + — ► pn data using Eqs. (46) 
and (47). The dashed lines correspond to the upper and lower bounds from the HyperCP result. 



therefore, the 7 8 contributions are important only in the kaon cases. We also note that nipping the 
signs of A pn o and B pn o in Eq. (47), whose overall sign is not fixed by experiment, would prevent the 
cancelation in the hyperon case from occurring and thus resulting in rates much above the bounds. 

It turns out that the imaginary part of the penguin amplitude is sufficient to eliminate these 
scalar particles as candidates for the HyperCP events, as it cannot be canceled by the 4-quark 
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amplitudes [29], having a size of 

\Img H \ ~ 5.8 x 10" 7 , (50) 

much larger than allowed by Eq. (9) with ImM^ = 0. The scaling of the penguin amplitude to 
the 5-meson system is also incompatible with the b — > sX bound. In the SM 

9'n = E m XsV qb ~ -1-7 x 10- 4 , (51) 

V q=u,c,t 

which is much larger than allowed by Eq. (20). In the 2HDM, the relative size is also too large: 
Wh/9h\ ~ \ v tb/V td \ 4 

Both of these problems are associated with a structure in which the Higgs-penguin amplitude 
is dominated by diagrams with up-type quarks and W bosons in the loops. It may be possible to 
remedy these problems in models with additional contributions to the penguin, for example, from 
super symmetric (SUSY) partners. If the penguin can be sufficiently suppressed, Eqs. (44) and (48) 
suggest that models in which kw ~ kc could satisfy the kaon bounds while being able to account 
for the HyperCP result. 



IV. PSEUDOSCALAR HIGGS BOSON 



We now consider the possibility that the new particle is a light CP-odd pseudoscalar, A, in the 
two-Higgs-doublet model. Specifically, we do so in types I and II of the model. 



A. Two-quark |A5| = 1 interactions 



The 2-quark flavor-changing couplings of A in the 2HDM are induced at one loop and have been 
evaluated in Refs. [31, 32]. The effective Lagrangian is the same in types I and II of the model and 
can be written as C As d in Eq. (lb), where 

9a = E mlV* d V qs ( + , (52 ) 

16^2 n 2 ^ q q q \ tan tan 3 /?/ v ; 

v q=u,c,t \ / 

with Ai^iq) being functions of m , m w , and m H +, whose expressions can be found in Ref. [31]. 
The leading-order chiral realization of Cj,sd is then 

C A = b D (B {h A , B}) + b F (B [h A , B])+b (h A ) (BB) + \fB, (h A ) + H.c. , (53) 



One could arrive at a similar conclusion about 7i in the 2HDM II by analyzing the decay tj —> ir°H, whose 
amplitude depends on the 4-quark parameters k d —k u [27]. Thus, from the 90%-C.L. bound B(r\ — > 7r°7i) < 
5 x 10~ 6 [30], one extracts \k d — k u \ < 0.45 for m H = 214.3MeV, which is incompatible with the limit derived 
from Eq. (43) plus the LEP constraint sin 2 (/3 - a) < 0.06 [21], namely \k d - k u \ = |2 cos(a - /?)/ sin(2/3)| > 1.9. 
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where 

h A = -2ig A (?hM? - tMh£) £ . (54) 

The leading-order diagrams for K — > it A and S + — > arising from £.4, plus £ s , are similar 
to those in the case of standard- model Higgs boson, displayed in Fig. I. The resulting amplitudes 
are 

M 2 ,(K + - ^] = -V2 M 2 ,(K a »M) = i 9j) m * ~ , (55) 
.M 2 ,(£+ - M) = *U ?S + - «U (D - F) - J* p T5 E+ . (56) 



B. Four-quark |A5| = 1 interactions 



The diagonal couplings of A to light quarks in the 2HDM are described by 

%A ~ %A 

^A qq = -qM-f 5 q— = -q L Mq R — + H.c. , 



where 



with 



q = ( u d s) T , M — diag(/ M m, l d m, l d m s ) , 



L = ~h = -cot (3 
I — — cot P , L = — tan /5 



in the 2HDM I 



in the 2HDM II . 



(57) 



(58) 



(59) 
(60) 



Since the Lagrangian for the quark masses is C q = —q L Mq R + H.c, the effect of C Aqq on interactions 
described by £ SjW can be taken into account using £ SiW and substituting M with MiA/v [10]. The 
resulting Lagrangians are 

Cf = (b D (B{M_,B}) + b F (B[M_,B]) + b (M_)(BB) + |/ 2 5 (M_>) ^ , (61) 



iA 



where 



Ci = 2rr s fB (h£M_?)— + H.c. , 



M_ = CMC ■ 



(62) 



(63) 



In addition, if the SU(3) singlet % is included in d$ by replacing E with £ exp(«-^/2/3r/ 1 //), the 
coupling of „4 to two gluons via the axial anomaly gives rise to [10] 



n \ ( 2 2 2 1 2^ 



(64) 
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which modifies the r\ x -A mixing generated by Cf. 
From dS 

, we derive the leading-order diagrams shown in Fig. 7 for K — > 7T.A, where 

rj = r) 8 cos 9 — 7] 1 sin , rj = rj 8 sin 6 + i] 1 cos . (65) 
The resulting amplitudes are 

+ Hs [{ 2m K + m l~ 3m 2 A ) c g -V8 (m 2 K - m 2 n ) s e ] 

[4l d m\ - (3l d + y ml] c + V 2 [2l d m\ + l u m 2 n - (l d + 2l u ) m 2 ] s 



x 



6 (rrii — m A ) v 
+ ils [{ 2m K + m l~ 3™>a) s g + V8 (m 2 K - ml) cj 

~ (3i d + Zj ml] s e -V2 [2l d m 2 K + l u m\ - (l d + 2l u ) m 2 ] c 



6(m 2 , - m A ) v 
h 8 ft, ~ Q i 2m K -ml- m 2 A ) ml 



where 



(66a) 



(m A — ml) v 
+ Hs [{ 2m K + m l~ 3m 2 A ) c g -V8 [m\ - m\) s e ] 
x [4/ d m\ - (31 d + l u ) mj] c + V2 [2l d m\ + l u mj - (l d + 2jJ m 2 ] s 

< o\[2 (m A — mi) v 
+ h 8 [{ 2m K + m l~ 3m A) s e + V8 (m 2 K - ml) c e ] 
x [4Z d m\ - (3l d + y ml] s e -^2 [21 d m\ + l u ml - (/, + 2/J mjj] c e 

6^ (m 2 A - m 2 ,) v 



cos 9 , s e = sin# , = — |m|- — . (67) 



The 7 8 contributions to this amplitude cancel completely, as already noted in Ref . [10] . Numerically, 
m ~ 819 MeV from fitting to the rj' mass after diagonalizing the r] 8 1 masses derived from the 
Lagrangians in Eqs. (28) and (64), and consequently 9 ~ —19.7°. 



A 



A 



A 



K ■ 



K - 



7T°, T], Vj 



K- 



K- 



FIG. 7: Diagrams contributing to K — > ttA arising from the 4-quark operators. The dots come from C 
in Eqs. (28) and (61), whereas the square vertices are from cffl in Eqs. (36) and (62). 



(A) 
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The leading-order 4-quark contributions to S + — > yx4 arise from the diagrams in Fig. 8 and can 
be expressed as 

M Aq ^ + - P A) = ip (A pA - B pAlh ) £+ , (68) 

where 

/ Kfl ( l d - Q ml 



2{m A — m^) f 

/ V { [Al d m 2 K - (3l d + g ml] c 2 e + ^2 [2l d m\ + l u ml - (l d + 2Z tt ) m 2 ] c^} 



_l_ - pi- 



2(m^ — m^) f 



/V { [^m 2 * - (3l d + g m 2 ] s 2 , - V2 [2Z d m| + / u m 2 - (Z d + 2/J m 2 ] Ce% } 

+ ^r~2 r\ ' ( 69 ) 



B 



fB pn o (l d ~l u ) 
pA 2(m A -ml)v 

f B o { [Al d m\ - (3l d + g ml] c 2 + ^2 [21 d m\ + Z u m 2 - (l d + 2l u ) m 2 ] c e s e } 



2 (m 2 — m^) v 



f B P ,° { [U d m 2 K - (3^ + g ml] sl-V2 [2l d m\ + l u ml - (l d + 2l u ) m 2 ] c^} 

o / 2 2 \ ' I' U J 

where y4 p7r o and -B p7r o are given in Eq. (46). We note that contributions with 7 8 or 7 8 appear only 
at next-to-leading order. 



C. Total contributions 



The total amplitudes for K — > ttA result from adding the contributions in Eqs. (55) and (66). 
If the CP-violating terms in the amplitudes are ignored, it is possible for the 2-quark and 4-quark 
contributions to cancel. We show this possibility in Fig. 9, where we plot the resulting branching 
ratios as functions of the charged-Higgs-boson mass for m A = 214.3 MeV and different tan f3 values 
in the 2 versions of the 2HDM. The total amplitude for S + — > pA is the sum of the contributions 



A 



A 



A 



7T ,77,?/ 



7T°, 77,77' 



71-°, 77,?/ 



FIG. 8: Diagrams contributing to S + — ► pA arising from the 4-quark operators. The square vertices come 
from C w in Eq. (36), whereas the dots are from the Lagrangians in Eqs. (28), (61), and (64). 
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2HDM I with tan/3 = 4 2HDM II with tan/3 = 0.9 
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ra#+ (GeV) m#+ (GeV) 

FIG. 9: Contributions of real parts of total amplitudes for K + — > ir + A (solid curve) and K$ — > tt°A 
(dashed curve) in the 2HDM to their branching ratios as functions of charged-Higgs-boson mass for = 
214.3 MeV and tan /3 = 4 (0.9) in type I (II) of the model. The horizontal lines indicate the upper bounds 
in Eqs. (7) and (8). 




225 250 275 300 325 350 375 400 



m H + (GeV) 

FIG. 10: Contribution of real part of total amplitude for S + — > pA to its branching ratio in the 2HDM I 
(solid curve) and II (dotted curve) as function of charged-Higgs-boson mass for = 214.3 MeV and 
tan (5 = 4 (0.9) in type I (II). The dashed lines indicate the bounds from the HyperCP result. 

in Eqs. (56) and (68). If the experimental values of A pn o and B p7T o in Eq. (47), as well as the middle 
value D — F = 0.25, are used in the total amplitude, the resulting branching ratios in the 2HDM 
are displayed in Fig. 10. 

We have found that only one of the kaon bounds can be satisfied if the HyperCP result is 
assumed to be mediated by A in the 2HDM. However, for certain tan/? and values near 

the ones indicated in Figs. 9 and 10, all the kaon and hyperon constraints can be nearly satisfied 
simultaneously. Part of the difficulty in satisfying all of the constraints lies with the vanishing 
of the K + and K s rates occurring at different m H + values, which is due to A4 iq (K s — > n°A) and 
— ReA44 q (K + — ► n + A) being unequal with l u ^ l d in Eqs. (59) and (60). We note that the situation 
is not much different if the signs of A pn o and B p7T o in Eq. (47) are both flipped. 

To summarize this section, we have found that it is possible for the real part of the penguin 
amplitude to cancel against the 4-quark amplitude to approximately satisfy the kaon bounds while 
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explaining the HyperCP observation with a 2HDM pseudoscalar. Unlike the scalar case, there is 
no free hadronic parameter at leading order in x?T in this case. The cancelation must happen as 
a function of the short- distance parameters that determine the size of the amplitudes. 

A feature shared by scalars and pseudoscalars in the 2HDM is that the imaginary part of the 
penguin amplitude is incompatible with the kaon bounds in Eq. (9) and has no counterpart that 
could cancel it in the 4-quark amplitude. A related problem is that the scaling of the penguin 
amplitude to the B system is also incompatible with observation. 

In view of these flaws, it is tempting to search for a model in which the penguin amplitudes are 
completely suppressed, and the 2HDM II seems to allow us to do that. In the 2HDM II the penguin 
amplitudes are proportional to l u , whereas the 4-quark amplitudes receive contributions from both 
l u and Id in Eq. (60). Thus the model in the large-tan /3 limit has l u — > 0. Unfortunately, in this 
limit Id induces 4-quark amplitudes resulting in 

£(£+ -> pA) 



B{K+ -> ir+A) 



0.025 , (71) 



which is inconsistent with Eqs. (9) and (13). In the 2HDM I, which has l u and Id given in Eq. (59), 
the 4-quark amplitudes alone yield 



B 4q (K+^n+A) 



0.53 (72) 



for all values of tan/3, which is consistent with Eqs. (9) and (13). However, in this case it is the 
penguin amplitude that eliminates the pseudoscalar as a possible HyperCP candidate. 

These results suggest the ingredients of a model that can satisfy all constraints. It is necessary 
for the penguin amplitudes to be dominated by additional particles, such as SUSY partners, in such 
a way that is not proportional to top-quark CKM angles. We have sketched a scenario where 
this happens in Ref. [33]. 



V. SUMMARY AND CONCLUSIONS 



We have summarized the existing constraints on the production of a light Higgs boson in kaon 
and 5-meson decays, as well as the implication of attributing the HyperCP events to the production 
of a light Higgs boson in hyperon decay. 

Production rates for such a particle in kaon and hyperon decays receive contributions from two- 
and four-quark operators that can be comparable in some cases. We have investigated the interplay 
of both production mechanisms with the aid of leading-order chiral perturbation theory. To this ef- 
fect, we have implemented the low-energy theorems governing the couplings of light (pseudo) scalars 
to hadrons at leading order in baryon %PT, generalizing existing studies for kaon decay. 

We first discussed the case of a scalar Higgs boson. We found that the leading-order amplitudes 
in both kaon and hyperon decays depend on an unknown low-energy constant 7 g , as well as known 
constants from the hyperon sector. This constant is connected to a weak-mass term in the chiral 
Lagrangian that can be rotated away for processes that involve only pseudo-Goldstone bosons and 
is, therefore, unknown. We applied our results to the process S + — > pX relevant to the HyperCP 
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observation of S + — > We showed that the two-quark contributions in the SM and its 

2HDM extensions are too large to explain the HyperCP observation. However, we also showed 
that there can be cancelations between the CP-conserving two- and four-quark contributions to 
this process that lead to a rate comparable in size to the HyperCP observation for both the SM 
and the 2HDM. Such cancelations occur for a certain range of known constants from the hyperon 
sector, the effect of 7 8 being small. In both cases, however, the two-quark penguin contribution 
has an imaginary (CP violating) part that is too large to be compatible with the HyperCP result. 
In the SM and in the 2HDM, the four-quark contributions have a CP-violating part that is much 
smaller than that of the penguin amplitude and hence these models are ruled out as explanations 
for the HyperCP observation. More general models with additional CP-violating phases may be 
able to address this issue. In addition, in these models the scaling of the two-quark operator to the 
B system is incompatible with the nonobservation of a light scalar in B decay. 

We then discussed the case of a pseudoscalar Higgs boson in the 2HDM. In this case we computed 
the leading-order amplitudes in %PT and included, as well, certain higher-order terms mediated by 
the r/ state. The resulting amplitudes for both kaon and hyperon decays do not depend on any 
unknown hadronic parameters. In particular, they do not depend on 7 8 , as observed in Ref. [10]. 
We then applied our results to the S + — > pA process. Once again we found that the real part of 
the amplitude can be consistent with the HyperCP observation for a certain range of parameters 
in the 2HDM (tan/5 and m H +), but that the imaginary part of the penguin amplitude is too large. 
The scaling of the two-quark operator to the B system also produces aB^ X S A rate that is too 
large. Both of these problems can be solved in more general models that modify the phase and 
scaling with CKM angles of the two-quark operator. 

In conclusion, we have shown that it is possible to interpret the HyperCP observation as evidence 
for a light Higgs boson, although it is not easy to arrange this in a model. Typical Higgs-penguin 
operators have three problems: 

(a) if they have the right size to fit the HyperCP observation, they induce K — > nX at rates 
larger than the existing bounds; 

(b) if they are dominated by loop diagrams involving up-type quarks and W bosons, they have 
a CP phase that is too large; 

(c) if they are dominated by loop diagrams involving up-type quarks and W bosons, their scaling 
to the B system is incompatible with the nonobservation of B — > X S X. 

We have found in this paper that (a) can be solved in some cases by the addition of the effects of 
four-quark operators. We have suggested that more general models may be constructed to solve (b) 
and (c). To show that this is possible, we have constructed a specific example in Ref. [33]. 

Disregarding existing bounds from kaon and P-meson decays, we have shown that many light 
Higgs bosons have couplings of the right size to explain the HyperCP observation. We think this is 
sufficiently intriguing to warrant a revisiting of the kaon and B decay results. In particular, the B 
factories are still operational and could reanalyze the very low invariant-mass region in their 
measurements of B — > X s fi + fi~ modes. The NA48 experiment might also be able to revisit the 
kaon modes. 
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